Report for Model 8

— for both jointly and separately

2024-08-19

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.07
## LDevsum       1.45       2.35
## dh0           1.09       1.31
## dh1           1.09       1.32
## dl0           7.54      16.72
## dl1           1.30       1.93
## dl2           6.79      15.95
## dl3           3.69       7.88
## 
## Multivariate psrf
## 
## 4.97

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1216.28632 22069.538
DIC3 1156.95606 22120.978
PWAIC 42.68576 234.347
WAIC 1183.68816 22141.168

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum      0.999       1.00
## LDevsum      1.001       1.01
## dh0          1.022       1.10
## dh1          1.028       1.12
## dl0          1.651       3.04
## dl1          1.401       2.27
## dl2          1.257       1.95
## dl3          1.446       2.33
## 
## Multivariate psrf
## 
## 1.47

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1213.74695 22043.9831
DIC3 1157.49543 22121.3879
PWAIC 43.02139 235.2914
WAIC 1184.60205 22142.1846